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We describe effects of nonlinearity on propagation of surface plasmon 
polaritons (SPPs) at an interface between a metal and an amplifying medium 
of the externally pumped two-level atoms. Using Maxwell equations we derive 
the nonlinear dispersion law and demonstrate that, the nonlinear saturation 
of the linear gain leads to formation of stationary SPP modes with the 
intensities independent from the propagation distance. Transition to the 
regime of stationary propagation is similar to the threshold crossover in lasers 
and leads to narrowing of the SPP spectrum. © 2009 Optical Society of 
America 

Delivery of electromagnetic energy from macro- to nanoscales and its direct generation 
in nanostructures are the challenging problems of nanophotonics. They are particularly im- 
portant for plasmonics and metamaterials, which promise subwavelength localization and 
advanced control of light in photonic nanocircuits. However, metal structures typically suffer 
from large intrinsic (Ohmic) losses hampering attractive applications. On the other hand, 
the electromagnetic field confinement and enhancement associated with surface plasmons 
have been proposed to control excitation of active emitters and achieve a plasmonic laser 
(spaser) [1,2]. Propagating surface plasmon polaritons (SPPs) at an interface with an ampli- 
fying dielectric material have been studied theoretically, see, e.g, [3,4] and references therein, 
and practically demonstrated using optically pumped dyes [5-7] and erbium doped glass [8]. 
In particular, a distinct threshold in the dependence of the output SPP intensity on gain 
and the simultaneous significant narrowing of SPP spectra have been reported inferring 
signatures of stimulated SPP emission. [7]. 
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Theory of SPPs interacting with an amplifying medium has been so far limited to the linear 
case [1,3,4,6,7]. This approximation describes the SPP modes, where the linear growth of the 
amplified SPPs is not balanced by the nonlinear losses. In the linear theory plasmons with 
the real propagation constant (in the SPP case) or the real frequency of the oscillations (in 
the spaser case) can be found only at the threshold, where the linear loss is exactly balanced 
by the linear gain [3]. However, in practise, the gain-loss balance should be maintained even 
above the threshold providing that nonlinearity and hence the SPP intensity are accounted 
for. This is because the amplification of SPPs is expected to be fully compensated by the 
nonlinear losses ensuring existence of SPPs with the real propagation constants. 

Our main aim here is to present a theory of SPP amplification accounting for saturation 
of the linear gain by the nonlinear losses. This theory predicts appearance of SPP modes 
with the stationary values of the amplitudes reached above the threshold and reveals crucial 
role of nonlinearity in shaping spectra of the excited SPPs. We start our analysis from the 
time independent Maxwell equations for TM waves 

8 ZZ E' X - d zx E' z = -D x , (1) 
d zx E' x - d xx E' z = D z , (2) 

where c is the vacuum speed of light, z is the coordinate along the interface and the x 
direction is orthogonal to it, both are measured in the units of 1/k = \ vac /(2ir), where 
\ vac is the vacuum wavelength. For the constitutive relation we assume D — (e + 7|-E ,/ | 2 )£' / , 
where e = e' + ie" is the linear permittivity and 7 = 7' + 27" is the nonlinear susceptibility. 
Below we use e and 7 with subscripts d and m when referring to the dielectric (x > 0) 
and the metal (x < 0), respectively. The metal is assumed linear, 7 m = 0. Amplification in 
the dielectric is described using the two-level model with relatively small SPP intensities. 
Susceptibility of the two-level atoms is x(^) = 01 (5 — i)/ (r 2 + 5 2 ) [9], where 5 is the detuning 
of the SPP frequency from the atomic resonance frequency and T = (1 + \E'\ 2 ^E^ 2 ) 1 ' 2 is the 
dimensionless intensity dependent linewidth. 5 and T are both normalized to the physical 
transition linewidth. a is the dimensionless gain per unit length. We assume below that 
the SPP intensity is less than the saturation intensity \E*\ 2 , a particular value of the latter 
depends on a material choice and is not important here. The effect of the nonlinear saturation 
of the linear gain discussed below should not be confused with and attributed to \E*\ 2 . 

We seek solutions of Eqs. (1-2) in the form E'(x,z) = E*E(x) exp(i(3z), where (3 is the 
propagation constant. Expanding T and \ into the Taylor series in \E\ 2 we find = e?, + 
a(5 + i)/(l + S 2 ), 7,2 = a(i — 5) / (1 + 5 2 ) 2 , where is the dielectric constant of the background 
material hosting the two-level atoms. For a > we have linear gain and nonlinear absorption, 
which both are maximal at the line center 5 = 0. We are not taking into account the metal 
dispersion, assuming that it is negligible relative to the dispersion introduced by the strong 
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two-level resonance. Linear limit of the theory developed below admits arbitrary complex 
/5's. However nonlinear results require the stationarity of the SPP intensities with respect to 
the propagation coordinate z, which is achieved due to balance between all the loss and gain 
mechanisms. Formally the balance condition is expressed as Im/3 = 0. 

The exponentially decaying for x < solutions can be readily found, since the problem is 
linear inside the metal: E x = B e q,nX , E z = (iq m //3)B e qmX , Re(q m ) > 0. Here q 2 m = /3 2 — e m 
and \B\ 2 is the SPP intensity on the metal side of the interface. The system of equations we 
need to solve for x > is 

(3 2 E X + i(3d x E z = [e d + ld {\E X \ 2 + \E Z \ 2 )} E x , (3) 
i/3d x E x - d xx E z = [e d + ld [\E X \ 2 + \E Z \ 2 )\ E z , (4) 

The boundary conditions require continuity of E z and D x at x — 0. Assuming that E x q, E z q 
are the field components on the dielectric side of the interface and knowing solutions inside 
the metal we express the boundary conditions using only the fields in the dielectric 

pe m E z0 = iq m [e d + ^ d (\E xQ \ 2 + \E z0 \ 2 )] E x0 , (5) 

Solving Eqs. (3,4) perturbatively under the assumption that nonlinear terms 7d|£'| 2 are small 
we found 

E x = Ae~^ x {l + w xld \A\ 2 e~ 2xR ^ + 0(\ ld \ 2 )}, (6) 
E z = %Ae-™*{\ + Wyld \A\ 2 e- 2xRe ^ + 0(| 7d | 2 )}. (7) 

Here q 2 d = f3 2 — e d , Re(q d ) > 0, and w x>y are some constants not shown here. \A\ 2 characterizes 
the SPP intensity on the dielectric side of the interface. Substituting the above solutions into 
Eq. (5), we find the dispersion law for SPPs accounting for nonlinearity, losses in metal and 
gain in the dielectric 

e m q d + e d q m =~f d \A\ 2 F + 0(\-f d \ 2 ). (8) 
The constant F is given by 



F.^ + ljx (9) 
+ q m e d + 2Re(q d )(P 2 e m /e d - q m q d ) 



ARe(q d )(Re(q d ) + q d ) 

Dispersion of nonlinear SPPs in the absence of gain and loss has been previously derived, 
e.g., in Ref. [10]. 

For \A\ =0, Eq. (8) transforms into a well known linear dispersion law for SPPs, which is 
readily resolved with respect to (3: (3 = (3i = , / - e ^_ € e d . Practically, SPPs can be excited either 
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directly by the two-level emitters (dye molecules, quantum dots, etc) or externally coupled 
into the system, e.g., through a prism or a grating. Depending on the spatial variations of 
the emitter density and on the excitation type (optical, electric or chemical pumping), the 
population inversion and hence the gain coefficient a may vary with the distance from the 
interface, see, e.g. [4]. We will not consider these effects in order to focus on the role of 
nonlinearity and to derive closed analytical expression for the nonlinear SPP dispersion. To 
obtain physical estimates for our dimensionless calculations, we used e& = 1.8 (polymer) and 
e m = — 15 + i0.4 (silver at A mc = 530 nm). For these parameters, without the resonant atoms 
(a = 0) the characteristic SPP propagation distance is (klmPi)" 1 ~ 30/im. 

If we neglect the effect of nonlinearity, the condition Imfy = corresponds to the lossless 
SPP propagation and it is achieved at the gain threshold a = ao: a o — Yrr[{\ e m\ 2 — 2efie&5) — 
{(l e m| 2 — 2e^e b 5) — 4(e^) 2 (e b ) 2 (1 + 5 2 )} 1 / 2 ]. As expected, the lowest gain a = a min required 
for the lossless propagation happens at the exact resonance (5 = 0): a min ~ 0.00575. For 
example for a = 1.5a m j n and 2a miri the characteristic SPP gain length {klmf3i)~ l ~ 65/im 
and 30/zm, respectively. The gain length is defined as the propagation distance at which the 
SPP amplitude increases by a factor e. Intensity of the nonstationary (Im/3 ^ 0) SPPs in 
the linear case can be easily calculated: // ~ e -2z/mft rj-^g ^ypj ca j dependence of l\ on 5 is 
shown in Fig. 1(a) (line 2). Lines 2,3 in Fig. 1(b) show how the full width half maximum 
(FWHM) of I i (5) varies with the gain parameter a for two different propagation distances. 
One can see that the spectrum quickly narrows as the gain is increased, but kept below the 
threshold (a < a m i n ). Close to and above the threshold the narrowing continues, but at a 
much slower pace. With the increase of the propagation distance and for the fixed gain, the 
spectrum also narrows (cf. lines 2 and 3 in Fig. 1(b)), since the spectral components of SPP 
modes near the line center are stronger amplified and hence become dominant. 

The influence of nonlinear effects on the SPP propagation constant can be derived by 
solving Eq. (9) with respect to (3 and demanding Im/3 = 0. We assume, that the right-hand 
side of Eq. (9) can be treated as a perturbation and find: 

P = P l + MA\\^ ^ff\ + f\ (10) 

2e 2 d Re(q d ) + q d 

where qd, m inside (3 n i are calculated for j3 = j3\. Above the threshold (a > a m i n ), the linear 
growth of the SPP intensity is saturated by the nonlinear absorption and as a result the 
intensity quickly attains a stationary value. To calculate the spectral and other characteristics 
of the stationary (saturated) SPPs, we use Eq. (10) and impose the condition Imj3 = 0. 

In the linear approximation, Imfti = implies a = ao (see above) and the corresponding 
real propagation constant is 0i(ao) = Po- Expanding (3 in Eq. (10) into the Taylor series in 
(a — ao), i.e. close to the threshold for a given 5, we find (3 = /3q + (a — ao)d a (3i + /3 n ;|A| 2 . 
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Then Imf3 = gives the intensity I s = \A S \ 2 of the stationary SPPs: 

I s = (a- a )lmd a Pi/(-lm/3 n i), (11) 

where d a /3i and [3 n i are calculated for a = ao- Fig. 2 shows dependencies of I s on gain for 
several values of 5. Naturally, one can see that for higher gain, SPPs within the wider fre- 
quency interval around the resonance cross the threshold, leading to the gradual broadening 
of the spectra of the stationary SPPs (see Line 1 in Fig. 1(b)). 

Fig. 1(a) compares the spectral intensity profiles of the linear SPPs below threshold and of 
the nonlinear saturated ones above the thresholds. Also, Fig. 1(b) compares the linewidth of 
these two SPP families. One can see that the spectra of stationary SPPs above the threshold 
are much narrower for small deviations from a m i n , than spectra of the linear SPPs. For prop- 
agation distances of few gain lengths (klm/3i)~ l (e.g., 100/im) we shall expect that the SPPs 
should achieve their stationary saturated intensities. Therefore, the threshold crossover, if 
observed at sufficiently long distances, should be accompanied by a marked spectral nar- 
rowing, which has indeed been reported in the recent experiments [7]. Note the obvious 
qualitative agreement between our Figs. 1(a) and 2 and the experimental Figs. 2(a) and 2(b) 
in Ref. [7]. 

For the propagation constant of the stationary saturated SPPs we thus have /3 S = (3q + 
(a — a )Red a /3[ + I s Rej3 n i. Fig. 3 compares real and imaginary parts of the above (3 S with 
Pi. It shows that there exists a marked difference of the dependencies of the propagation 
constant from the gain parameter in the linear and nonlinear cases. We have independently 
cross-checked the SPP profiles and propagation constants using numerical shooting method 
applied directly to the Maxwell equations. Good agreement between analytical and numerical 
results for a < 2a m i n made it unnecessary to present the numerical results in the context of 
this paper. 

In summary, we have presented a theory of SPP amplification in the presence of nonlinear 
gain saturation as it happens in the real-world systems. We have calculated saturated inten- 
sities of SPPs above the amplification threshold and quantitatively described transition from 
the broad Lorentzian spectra of the nonstationary amplified SPPs to the narrow spectra of 
the stationary saturated SPPs above the threshold. Our analytical results are expressed in 
terms of complex dielectric permittivities of the materials involved and can be readily used 
in design of plasmonic devices. 

The work of A. Zayats has been supported in part by EPSRC (UK) and by EC FP6 
project PLASMOCOM. The authors are grateful to A. Krasavin for discussions. 
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Fig. 1. (Color online) (a) Normalized SPP intensity as a function of detuning 6. Line 1 is the 
spectrum of stationary (saturated) SPP, I s (6): a = l.la m j n . Line 2 is the spectrum of the 
amplified linear SPP, Ii(S), below the threshold (propagation distance 60/im, a = 0.8a m i n ). 
(b) FWHM of the SPP spectra vs gain. Line 1 is for the stationary saturated SPPs, I s . Lines 2 
(propagation distance 60/im) and 3 (distance 100/xm) are for the nonstationary linear SPPs, 



Fig. 2. (Color online) The dependence of the intensity of the stationary SPPs above the 
threshold on the gain parameter a/a m i n . Line 1 corresponds to 6 = 0, Line 2 to 6 = 0.4 and 
Line 3 to 5 = 0.8. 



Fig. 3. (Color online) The dependence of the real (a) and imaginary (b) parts of the SPP 
propagation constants on the gain parameter for 6 = 0. Solid red lines correspond to the 
stationary nonlinear SPPs (j3 s ) and the dashed black lines to the linear SPPs ($). 
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